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ABSTRACT 

The  power  spectrum  of  the  horizontal  component  of  the  solar  convective  velocity  field  has  recently  been 
estimated  from  observations  of  the  Doppler  shifts  of  surface  flows  at  and  near  disk  center  (Chou  et  al.).  From 
their  analysis,  the  authors  assert  that  “there  is  no  evidence  of  apparent  energy  excess  at  the  scale  of  meso- 
granulation.”  We  show  in  this  paper  that  their  conclusion  is  incorrect  and  that  the  shape  of  the  observational 
spectrum  does  indeed  confirm  the  presence  of  both  supergranules  and  mesogranules  in  the  solar  convective 
flow. 

To  establish  this  claim  we  have  extended  existing  kinematic  models  of  convection  at  the  solar  surface 
(Simon  and  Weiss;  Simon  et  al.)  and  have  introduced  power  spectra  diagnostics.  We  find  that  models  with 
supergranule  cells  alone  do  not  produce  spectra  that  match  the  observations,  but  if  mesogranules  arc  included 
then  there  is  excellent  agreement  between  the  model  and  observational  spectra  when  the  model  parameters  are 
chosen  to  be  consistent  with  proper  motion  and  Doppler  measurements.  Since  the  magnitudes  of  our  model 
spectra  are  consistent  with  earlier  Doppler  measurements  of  the  rms  horizontal  velocity  field,  but  are  signifi¬ 
cantly  less  than  those  of  the  observed  spectrum,  we  must  call  into  question  the  analysis  of  Chou  et  al. 

Subject  headings:  convection  —  Sun:  granulation 


1.  INTRODUCTION 

With  the  development  of  local  correlation  techniques  for 
measuring  the  proper  motions  of  granules  (November  1986; 
November  &  Simon  1988),  it  became  possible  to  observe 
directly  horizontal  flows  at  the  solar  surface.  For  example,  the 
existence  of  mesogranular  convection,  originally  deduct  from 
Dopplergrams  of  the  vertical  velocity  (November  et  al.  1981), 
has  been  unambiguously  confirmed  by  observations  of  ~6 
Mm  structures  in  horizontal  flow  and  divergence  maps  derived 
from  proper  motion  measurements  (Title  et  al.  1986;  Simon  et 
iL  1988;  November  1989).  Subsequently  it  has  been  shown 
that  the  observed  surface  velocity  fields  resulting  from  super- 
granules,  mesogranules,  and  exploding  granules  can  be 
modeled  surprisingly  well  by  a  collection  of  irrotational  and 
axisymmetric  velocity  sources  of  varying  sizes  (Simon  &  Weiss 
1989;  Simon  et  al.  1991). 

Recently,  power  spectra  of  the  horizontal  and  vertical  solar 
surface  velocity  fields  have  been  estimated  by  analyzing  spectra 
of  the  line-of-sight  velocity  fields  in  a  series  of  256"  x  240" 
Doppler  images  displaced  sequentially  from  the  disk  center 
(Chou  et  al.  1991,  hereafter  Paper  I).  The  images  were  filtered 
in  the  frequency  domain  by  excluding  power  with  periods  less 
than  12.8  minutes  to  eliminate  p-mode  contamination.  Anom¬ 
alous  Doppler  line  shifts  in  regions  of  high  magnetic  field  were 
also  removed.  With  a  1"  image  resolution  the  resultant  power 
spectra  span  the  range  of  wave  numbers  from  k  —  0.042  rad 
Mm’‘  (/  ==  2n/k  ~  150  Mm)  to  k  ~  4.2  rad  Mm'*  {/.  ~  1.5 
Mm)  and  easily  include  the  scales  corresponding  to  super- 

'  Operated  by  the  National  Optical  Astronomy  Observatories  for  the 
Association  of  Universities  for  Research  in  Astronomy,  Inc.,  under  contract 
with  the  National  Science  Foundation.  Partial  support  for  the  National  Solar 
Observatory  is  provided  by  the  USAF  under  a  Memorandum  of  Understand- 
i»f  with  the  NSF. 


granules  and  mesogranules.  When  plotted  on  a  log-log  scale 
the  power  spectrum  of  the  observable  component  of  the  hori¬ 
zontal  velocity  field  (Fig.  7  of  Paper  1)  can  be  characterized  as 
generally  decreasing  in  magnitude  from  k  ->  0.1  rad  Mm' '  to 
k  ~  4.0  rad  Mm ' '  with  a  peak  in  the  range  k  ~  0. 1 35-0.304 
rad  Mm' '  (corresponding  to  the  supergranule  length  scale  of 
~  30  Mm)  and  an  average  slope  beyond  the  peak  of  ~  -  1.64. 
A  bend  in  the  spectrum  occurs  at  k  --  1.05  rad  Mm' ‘  (/.  ~  6 
Mm),  followed  by  a  steeper  decrease  with  slope  -  4.0.  Since 
they  found  no  peak  in  the  spectrum  corresponding  to  the 
mesogranular  length  scale  (^6  Mm),  the  authors  concluded 
that  “there  is  no  evidence  of  appprent  energy  excess  at  the 
scale  of  mesogranulation,"  thereby  contradicting  the  results 
from  earlier  Dopplergram  and  granule  proper  motion  observa¬ 
tions. 

In  this  paper  we  show  that  the  observed  power  spectra  pre¬ 
sented  in  Paper  I  are  consistent  with  a  distinct  meso- 
granulation  convection  scale  and,  in  fact,  are  just  what  we 
would  reasonably  expect  from  previously  obscrvwi  surface  flow 
structures.  We  do  this  by  extending  the  kinematic  models  of 
the  convectively  driven  surface  velocity  fields  (Simon  &  Weiss 
1989;  Simon  et  al.  199!)  to  include  many  supergranules  and 
mesogranules  covering  a  300"  x  300"  area  at  the  solar  disk 
center  (§  2).  To  emulate  closely  the  spatial  resolution  and  field 
of  view  of  Paper  I’s  observations  we  project  the  observable 
component  of  the  horizontal  velocity  field  from  the  center  of 
our  model  flow  field  onto  a  256"  x  256"  observational  plane 
with  r  resolution.  Computing  the  model  power  spectra  we 
demonstrate  good  agreement  with  the  shapes  of  the  observed 
spectra  of  Paper  I  over  ranges  of  the  model  parameters  that  are 
physically  reasonable  (§  3).  Comparing  magnitudes  of  the 
model  and  observational  spectra  we  find  some  nontrivial  dis¬ 
crepancies  which  we  discuss  in  §  4  along  with  our  conclusions. 
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‘  2.  MODEL  OF  CONVECTIVELY  DRIVEN  SURFACE  FLOWS 

The  surface  flow  model  is  composed  of  a  collection  of  inde¬ 
pendent  velocity  source  functions  distributed  in  an  x  -  y 
Cartesian  coordinate  system  with  each  source  function  rep¬ 
resenting  the  contribution  of  a  single  supergranule  or  meso- 
granule  convection  cell.  Simon  &  Weiss  (1989)  have  shown 
that  for  each  source,  an  irrotational  velocity  field  v  derived 
from  a  Gaussian  velocity  potential  <(>,  where  r  =  —  is  a 
good  approximation  to  the  observed  horizontal  velocity.  We 
therefore  choose  a  velocity  potential  of  the  form 

<P(x,  y)  =  X  y)  +  X  y)<  (•) 

t=i 

where  N,  is  the  number  of  supergranules  in  the  domain  and 
the  number  of  mesogranules.  The  supergranule  (s)  and  meso- 
granule  (m)  potential  functions  are  given  by 


y)  =  ^  exp  |- 
y)  =  exp  I 


-[(x-Xi)^-F(y  -  yj* 


-[(x-x/-^(y-y/ 


where  (V,  R)  determine  the  amplitude  and  spatial  size  of  a 
velocity  source,  and  (xj,  yj  or  (xj,  y^)  its  center.  Sizes  of  the  N, 
supcrgranules  and  mesogranules  are  chosen  to  satisfy 
normal  distributions  of  variances  dR,  and  d/t.  about  the 
"average  values  <R,>  and  <R«>,  respectively.  The  normal  dis¬ 
tributions  are  slightly  truncated  by  placing  lower  bounds  R^'" 
(supergranule)  and  RJ*"  (mesogranule)  on  the  source  sizes.  We 
do  not  include  exploding  granules  in  our  model  for  two 
reasons:  (1)  the  spectra  of  Paper  I  were  filtered  to  exclude  all 
power  with  variations  less  than  12.8  minutes,  and  this  exceeds 
the  lifetimes  ( ~  10  minutes)  of  most  exploders  (Title  et  al.  1989), 
and  (2)  the  velocity  field  of  a  long-lived  mesogranule  closely 


approximates  that  of  a  collection  of  short-lived  exploders  nor¬ 
mally  distributed  about  the  mesogranule  center. 

Ignoring  curvature  effects,  we  determine  the  position  of  the 
source  centers  on  a  300*  x  300'  (218  Mm  x  218  Mm)  rec¬ 
tangular  domain  fixed  to  the  solar  surface  in  the  following 
manner:  imagine  a  uniform  hexagonal  grid  covering  the 
domain  with  nearest  neighbor  spacing  d,.  A  supergranule 
source  is  placed  in  each  hexagonal  cell  with  the  center  of  the 
source  (x^,  yJ  uniformly  randomly  distributed  within  a  circle  of 
radius  p,  about  the  cell  center.  The  centers  of  the  N.  meso¬ 
granule  sources  (xj,  y^)  are  uniformly  randomly  distributed 
over  the  entire  domain  with  the  average  number  of  meso¬ 
granules  per  supergranulc  being  N„,  where  =  SJN,.  AO 
source  functions  of  a  given  variety  are  placed  such  that  any  two 
source  centers  are  separated  by  at  least  a  distance  equivalent  to 
their  combined  size.  [We  have  also  computed  a  normal  dis¬ 
tribution  of  mesogranules  about  each  supergranule  center. 
For  variances  of  order  {dJ2)^  the  results  of  our  analysis  do  not 
differ  significantly  from  those  results  obtained  assuming  a 
uniform  random  distribution.  Hence  we  consider  only  the 
latter  distribution  in  the  remainder  of  this  paper]. 

To  compare  our  model  predictions  with  the  observations  of 
Paper  I  we  compute  from  the  velocity  potential  (eq.  [1])  the 
obMrvable  component  (r^o)  of  the  horizontal  velocity  (ra),  Le, 
that  component  that  would  contribute  to  a  line-of-sight 
Doppler  velocity  measurement: 

|r-r,l  • 

where  r  =  x«,  +  ye,  and  r,  is  the  location  of  the  disk  center  as 
V  ewed  from  Earth,  coincident  with  the  domain  center  for  the 
models  presented  here.  We  then  project  Vno  onto  a  256’  x  256* 
observational  plane  with  1*  resolution  also  spanned  by  an  x-y 
Cartesian  coordinate  system. 

In  Figure  1  we  show  contour  plots  of  Vuo  in  the  observa¬ 
tional  plane  corresponding  both  to  an  “idealized”  regular 


r-N  ■'* 


^  fell' 
f 


Fio.  Id 


rnj.  la  rM.  lu 

Fxj.  t. — Contours  of  the  observable  component  of  the  horizontal  velocity  field  Sgg  derived  from  the  solar  surface  flow  model  in  a  256'  x  256’  observation*| 
plane  centered  at  the  solar  disk  center  showing  (a)  only  supergranules  in  a  regular  hexagonal  array  and  (6)  a  full  complement  of  supergranules  and  mesograo 
randomly  distributed  as  described  in  {  2.  Solid  (dotted)  Unes  represent  positive  (negative)  values.  See  text  for  values  of  the  model  parameters. 
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of  identical  supcrgranules  (Fig.  la)  and  the  full  "real” 
pgdel  with  all  scales  and  variations  present  (Fig.  16).  Figure  la 
b  generated  with  all  the  parameters  set  to  zero  except  for 
j  *  40"  (29  Mm),  V.  =  7.5  x  lO'*'  min'*  (0.91  km  s-‘X  and 
/gs  =  12"  (8.7  Mm).  Figure  16  is  generated  with  the  same 
r  and  <R.>  values  and  has  p,  =  15'  (11  Mm).  dR.  =  5'  (3.6 
gmi.  ^  y  (2.2  Mm).  N.  =  12.  F.  =  8.0  X  10’ *'  min'' 
«97  km  s'  ‘X  <R.>  =  3"  (12  MmX  dR.  =  2"  (1.5  MmX  and 
=  r  (0.73  MmX  Where  possible,  the  parameters  for  the 
jjP  model  have  been  chosen  to  be  consistent  with  values  deter- 
from  proper  motion  measurements.  In  particular, 
(R^,  <(R.)  are  close  to  values  found  in  the  measure- 

^eots  of  Simon  &  Weiss  (1989)  and  Simon  et  al.  (1991). 
further  observational  constraints  imposed  on  the  integrated 
power  of  the  spectrmn  will  be  discussed  in  §  4. 

3.  MODEL  POWER  SPECTRA 

power  spectra  of  the  observable  component  of  the  horizon- 
lul  velocity  are  constructed  using  the  standard  periodogram 
luetbod  (Press  et  al.  1986)  generalized  to  two  dimensions. 
Before  Fourier  transforming,  a  Hanning-type  window  function 
u  applied  to  the  spatial  field  in  a  25'  terder  region  at  the 
boundaries.  Labeling  the  two-dimensional  power  spectrum  as 
where  k,  and  k,  are  discrete  wave  vectors,  the  one- 
gimensional  power  spectrum  is  defined  as 

^  *,  t. 

where  k  =  mAk,  m  =  0, 1, 1 ....  128,  and  Ak  =  3.4  x  10"*  rad 
Mm'‘. 

We  consider  first  a  simple  model  comprised  only  of  super- 
granule  cells.  Figure  2  displays  three  model  power  spectra  cor¬ 
responding  to  different  values  of  the  size  distribution  variance 
dR,  and  maximum  cell  center  offset  />,  but  with  the  common 
parameter  values  d„  V„  <R,>,  and  Rf**  from  Figure  16.  It  is 
dear  that  (a)  if  the  supergranules  are  of  uniform  size  {dR,  =  0) 
and  fixed  to  the  centers  of  the  hexagonal  lattice  {p,  =  OX  the 
spectrum  is  peaked  in  the  range  k  0.22-0.29  rad  Mm~’.  (6) 
Randomizing  the  positions  of  the  supergranule  centers  (p,  = 
15")  about  the  centers  of  the  hexagonal  lattice  celk  broadens 
this  sharp  peak,  (c)  Further  broadening  occurs  when  the  sizes 
of  the  supergranules  are  distributed  normally  with  a  nonzero 
variance  (dR,  =  5")  resulting  in  a  spread  of  the  peak  power  to 
smaller  wavenumbers.  In  all  three  cases,  however,  the  spectrum 
drops  off  sharply  when  k  >  0.5  rad  Mm  ~ '  and  is  very  different 
from  the  observed  spectrum  of  Paper  I  (the  dot-dash  line 
superposed  on  Fig.  2)  which  is  seen  to  have  a  clear  shoulder  at 
k  1.0  rad  Mm~ '.  (Note  that  we  also  include  as  a  dotted  line 
the  spectrum  of  Paper  I  scaled  by  a  factor  0.13  so  that  it 
coincides  more  closely  with  the  model  spectra.) 

The  effect  of  adding  mesogranulc  velocity  sources  into  the 
model  is  illustrated  in  Figure  3.  Three  model  power  spectra  are 
shown  using  the  supergranule  parameters  and  common  meso- 
granule  parameters  Al.,  y„,  <R.>,  and  R;'*  from  Figure  16, 
but  with  different  values  of  the  size  distribution  variance  dR.. 
Examining  the  power  spectrum  for  mesogranules  that  are 
uniform  in  size  (dR.  =  0,  solid  tine),  we  see  that  the  sharp  decay 
for  k  >  0.5  rad  Mm ' '  that  characterized  the  supergranule- 
only  spectrum  is  replaced  by  a  gentle  mesogranular  peak 


Woveienqth  X  (orcsec) 

80  60  40  20  10  8  6  4  3 


Fmj.  2. — One-dimeasiona]  power  spectra  of  the  model  velocity  field  for 
different  values  of  the  supergranule  distribution  parameters  with  no  meso¬ 
granules  present  The  top  axis  shows  the  wavelen^  (il)  corresponding  to  the 
spectral  wavenumber  {k).  See  key  and  text  for  parameter  values.  For  compari¬ 
son  we  also  show  the  spectrum  of  Paper  I  as  published  {dol~dash  line)  and  when 
scaled  by  a  factor  0. 1 3  [dotted  biel. 


between  k~0.5  rad  Mm*'  and  k~  1.2  rad  Mm“‘.  If  we 
introduce  a  nonzero  variance  to  the  mesogranule  size  distribu¬ 
tion  (dR.  =  ITO,  short  dashed  line)  the  peak  flattens  into  a 
“shelf”  with  a  shoulder  at  k  —  1.0  rad  Mm*’.  Increasing  the 
variance  further  (dR.  =  2T0,  long  dashed  line)  steepens  the 
mesogranular  shelf  and  results  in  a  log-log  spectrum  that 
decreases,  on  average,  with  a  slope  — 1.4  from  the  super- 
granule  p^  to  the  shoulder  at  k  ~  1.0  rad  Mm*  ’.  Beyond  the 
shoulder  the  slope  steepens  to  4.3.  This  latter  spectrum 
resembles  closely  the  shape  of  the  observed  spectrum  of  Paper 
I  {dot-dash  line)  and  in  addition  matches  satisfactorily  the  mag¬ 
nitude  of  the  0. 1 3  scaled  Paper  I  spectrum  {dotted  line). 

All  the  power  spectra  displayed  in  Figures  2  and  3  have  been 
computed  from  a  specified  configuration  of  supergranule  and 
mesogranule  velocity  source  functions  with  identical  random 


Wovelenqth  X  (orcsec) 

80  60  40  20  10  8  6  4  3 


Wovenumber  k  (rod  Mm"') 

Fic.  3.— Ooe-dimensioDal  power  spectra  of  the  model  velocity  field  for 
a  fixed  set  of  supergranule  parameters  but  varying  values  of  the  mesogranule 
distribution  parameters.  See  key  and  text  for  parameter  values.  The  scaled  and 
unsealed  spectra  from  Paper  I  are  also  showiL 
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Wovelenoth  X  foresee) 

80  60  40  20  10  0  6  4  3 


Wavenumber  k  (rod  Mm  ') 


Fk.  4. — The  ivente  at  the  ooe-dimensional  power  spectra  for  the  modd 
velocity  6eld  (joW  fine)  over  25  distioct  coofiguratioos  of  the  velocity 
source  functions.  The  short-dashed  line  above  and  the  long-dashed  line  below 
the  average  represent  the  maxiinum  and  minimum  power,  respectively,  found 
at  each  k  value  in  the  25  spectra.  Sec  text  for  values  of  the  m<^l  parameters. 
Again  we  show  the  scaled  and  unsealed  spectra  from  Paper  L 


distributions  of  sizes  and  positions  (except,  of  course,  when  the 
variance  was  zero).  To  estimate  the  statistical  reliability  of  our 
model  spectra,  we  used  a  random  number  generator  to 
produce  2S  different  configurations  with  the  same  model 
parameters  used  for  Figure  Ifc.  From  the  power  spectra  for  all 
these  configurations  we  then  determined  the  maximum, 
minimum,  and  average  power  for  each  disaete  wavenumber 
bin  and  plotted  the  result  in  Figure  4.  As  before,  we  superpose 
the  unsc^ed  and  scaled  observational  spectra  from  Paper  I 
(dot-dash',  dotted  lines).  The  distance  between  the  maximum 
curve  (short  dashed  line)  and  minimum  curve  (long  dashed  line) 
is  characteristic  of  the  uncertainty  in  any  power  spectrum  com¬ 
puted  from  a  distinct  configuratioa  Properties  of  the  average 
power  curve  (solid  Hne)  are  only  slightly  different  fi-om  the 
distinct  configuration  shown  in  Figure  3.  An  essential  point  to 
note  is  that  all  model  spectra  in  this  parameter  regime  retain 
the  signature  of  the  mesogranule  convection  in  the  form  of  a 
shoulder  at  Ji  ~  1.0  rad  Mm  ~ 


4.  DISCUSSION 

Comparing  the  shape  of  the  model  spectra  (Figs.  3  and  4) 
with  the  observational  spectrum  of  Paper  I  we  find  qualitative 
agreement :  at  low  h  both  spectra  have  a  peak  corresponding  to 
the  supergranule  spacing  that  gives  way  to  a  decreasing  func¬ 
tion  with  relatively  shallow  slope  interrupted  by  a  bend  at 
k  ^  1.0  rad  Mm'*  followed  by  a  more  strongly  decreasing 
function  at  hi^ier  k.  Values  of  the  logarithmic  slopes  and  the 
locations  in  k  of  the  supergranule  peak  and  mesogranule  break 
are  in  reasonable  agreement  given  the  statistical  deviations  of 
the  model  and  the  experimental  uncertainties  of  the  observa¬ 
tions. 

However,  comparing  the  magnitude  of  the  model  spectrum 
with  the  observational  spectrum  of  Paper  I  we  find  a  consider¬ 
able  discrepancy.  For  each  k  the  power  in  the  observational 
spectrum  is  7-8  times  larger  than  that  of  the  model  spectrum 


(Fig.  3).  In  our  model,  the  parameters  have  not  been  chosen 
simply  to  fit  the  shape  of  the  observed  spectrum,  but  rather  to 
be  consistent  with  the  parameter  values  determined  from 
proper  motion  measurements.  An  additional  constraint  on  the 
model  is  the  estimate  of  the  magnitude  of  the  horizontal  veloc¬ 
ity  field  derived  from  spatially  averaged  Doppler  velocity 
observations  (Hart  1956;  Simon  &  Leighton  1964;  GiovanelU 
1980).  By  Parseval’s  theorem,  the  integrated  power  spectrum  is 
equal  to  the  mean  square  amplitude  where  the  average 
is  over  the  entire  spatial  domain.  Integrating  the  average 
power  spectrum  of  Figure  4  we  find  an  rms  observable  com¬ 
ponent  of  the  horizontal  velocity  *  0.23  km  s'*. 

Converting  this  to  the  peak  horizontal  velocity  Vg  by  multi¬ 
plying  by  a  factor  of  2  (Giovanelli  1980)  we  obtain  Vg  %  0.46 
km  s  ~  *,  a  value  consistent  with  the  0.3-0.5  km  s  “ '  determined 
from  spatially  averaged  Doppler  velocity  observations. 

We  therefore  believe  that  the  discrepancy  between  the  model 
and  observational  spectra  lies  in  the  data  analysis  performed  in 
Paper  1.  We  have  discussed  with  the  authors  of  Paper  I  their 
inconsistent  labeling  of  graphs  (Chou  1991),  confusion  of  the 
horizontal  velocity  with  observable  component  of  the  horizon¬ 
tal  velocity  (Chou  1991),  and  the  absence  of  any  explicit  refer¬ 
ences  to  the  applied  window  functions  (Chou  1991)  in  the 
power  spectrum  decomposition  analysis.  These  problems  make 
it  impossible  for  us  to  reproduce  the  power  spectrum  of  the 
observable  component  of  the  horizontal  velocity  field  (Fi^  6  of 
Paper  I)  from  the  power  spectra  of  the  Doppler  velocities  on 
different  patches  of  the  solar  surface  (Fig.  2  of  Paper  I).  We 
assume  that  the  pdtvcr  spectrum  shown  in  Figure  6  of  Paper  1 
is  as  labeled;  i.c.,  the  power  spectrum  of  the  velocity  function 
defined  in  equation  (A2)  of  Papier  I.  If  we  estimate  the  rms 
velocity  associated  with  the  supergranule  peak  by  taking  the 
square  root  of  the  sum  of  the  power  in  the  wavenumber  range 
k~  0.14-0.30  rad  Mm"‘  (Parseval's  theorem),  we  then  find 
<p|o> 0.35  km  s“*  corresponding  to  a  peak  velocity 
Vg  ~  0.7  km  s"  *.  This  number  is  already  larger  than  the  value 
determined  by  earlier  Eioppler  velocity  measurements.  The  dis- 
crepiancy  is  made  considerably  worse  by  the  realization  that  a 
true  compiarison  requires  the  integration  of  the  power  spiec- 
trum  to  wavenumbers  characteristic  of  the  resolution  of  the 
Doppler  measurements;  namely,  to  values  k  s  1.0  rad  Mm  *■ 

In  summary,  we  have  shown  in  this  paper  that  a  model 
compiosed  of  both  supiergranule  and  mesogranule  convective 
surface  flows  produces  a  pmwer  spiectrum  of  the  azimuthal 
compxrnent  of  the  horizontad  velocity  that  is  consistent  with  the 
shapie  of  the  observed  spectrum  reported  in  Paper  I.  We  believe 
that  discrepancies  in  the  relative  magnitudes  of  the  model  and 
observed  spiectra  arise  from  errors  in  the  analysis  of  Paper  I- 
We  caution  that  our  model  is  not  unique  in  being  able  to 
reproduce  the  spiecific  shapie  of  the  observed  spjectrum,  and  we 
do  not  claim  to  have  isolated  the  optimal  set  of  parameters 
most  characteristic  of  the  physical  processes  involved  in  solar 
convectioiL  Rather,  we  have  determined  a  set  of  parameters  for 
supiergranule  and  mesogranule  flows  that  is  consistent  with  (D 
spratially  averaged  Doppler  velocity  observations,  (2)  proper 
motion  measurements  of  the  granulation,  and  (3)  (with  the 
caveats  discussed  above)  the  power  spectra  of  Paper  1.  We 
conclude,  therefore,  that  the  sp^ra  of  Paper  I  do  indeed  sho* 
the  presettce  of  energy  at  the  scale  of  the  mesogranules,  and  we 
disagree  strongly  with  the  authors’  assertion  to  the  contrary,  j 

We  wish  to  thank  Steve  Keil,  Larry  November,  and  Nige* 
Weiss  for  helpful  discussions. 
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